Reconstruction of motional states of neutral atoms via MaxEnt principle 
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We present a scheme for a reconstruction of states of quantum systems from incomplete 
tomographic-like data. The proposed scheme is based on the Jaynes principle of Maximum En- 
tropy. We apply our algorithm for a reconstruction of motional quantum states of neutral atoms. 
As an example we analyze the experimental data obtained by the group of C. Salomon at the ENS 
in Paris and we reconstruct Wigner functions of motional quantum states of Cs atoms trapped in 
an optical lattice. 



I. INTRODUCTION 



A reconstruction of states of quantum systems from experimental data represents an important tool for a verification 
I ' of predictions of quantum theory. It also allows us to check the fidelity of quantum state preparation as well to study 
the fidelity of processing of information encoded in states of quantum systems. Complementary, discrepancies between 
j— ( ■ estimated (reconstructed) states based on the measured data and theoretical predictions can serve as an indicator of 
^ \ various noise sources which occur during the quantum information processing or in the measurement of quantum states. 
■ Without a priori assumptions about the character of physical processes and properties of reconstruction schemes the 
reconstruction cannot distinguish between imperfections related to an incoherent quantum state processing and non- 
ideal measurements Jlj] . Determination of limits for coherent control of quantum degrees of freedom or identification 
of sources of decoherence are essential for systems which are considered for quantum computing and information 
, processing ||. 

In atomic optics a highly coherent control of motional degrees of freedom has been achieved for trapped ions H 
and recently also for neutral atoms W>|. Cold atoms can be cooled into specific quantum states within micro- wells 
of an optical lattice which is induced by laser beams. Cold neutral atoms in optical lattices represent a promising 
system for quantum information processing. To verify a degree (fidelity) of coherent control over motional degrees of 

(— i . freedom of neutral atoms a reconstruction of their motional quantum states from measured data has to be considered. 
We develop a reconstruction procedure based on the Jaynes principle of Maximum Entropy {MaxEnt) (6[|| to achieve 
this goal. This scheme allows us to perform a state reconstruction from the experimental data obtained at the ENS 

^ '■ in Paris. 

In Section |l^ we present a brief description of the reconstructio n pro cedure of a density operator of a quantum system 



based on the Jaynes principle of Maximum Entropy. In Section III we utilize the MaxEnt principle for development 



of the reconstruction scheme of motional states of atoms and we perform numerical tests of our approach in Section 



. [II. Our reconstruction scheme is applied to the experimental data in Section [V 



II. MAXENT PRINCIPLE AND RECONSTRUCTION OF DENSITY OPERATORS 

Let us assume a set of observables G v (y — 1, . . . ,n) associated with the quantum system under consideration. 
This system is prepared in an unknown state p. Let us assume that from a measurement performed over the system 
meanvalues G v of the observables G v are found. The task is to determine (estimate) the unknown state of the quantum 
system based on the results of the measurement. Providing the set of the observables G v is not equal to the quorum 
(i.e., the complete set of system observables ||), then the measured meanvalues do not determine the state uniquely. 
Specifically, there is a large number of density operators which fulfill the conditions 

Trp {6} = 1, 

Tr (p {&} G v )=G v , v = l,2,...,n, (2.1) 

that is, the the normalization condition and the constraints imposed by the results of the measurement. To estimate 
the unknown density operator in the most reliable way we utilize the Jaynes principle of maximum entropy (MaxEnt 
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principle) |6|-p|,|l0| according to which among those operators which fulfill the constraints ( J2.l[ ) the most reliable 
reconstruction (estimation) p r is the one with the maximal value of the von Neumann entropy S(p) = — Tr(plnp): 

S(p r )= max S(p {d} )-yp {d} . (2.2) 



As shown by E.Jaynes || The operator which fulfills the constraints (2.1) and simultaneously maximizes the von 
Neumann entropy can be expressed in the generalized canonical form 



— exp(-5^A„G v ); (2.3) 



Z {G} 



Pr = 

where 

Z {d} (\ u ...,\ n ) = Tr[exp(-'E l KG v )], (2.4) 

is the generalized partition function and \ v are the Lagrange multipliers. The Lagrange multipliers A„ are chosen 



so that the density operator (2.3) fulfills the constraints (2.1) imposed by the results of the measurement. It is then 



obvious that the meanvalues G v which determine the density operator are related to the Lagrange multipliers via the 
derivatives of the partition function 

G v = Tr(p r G v ) = -^-lnZ {6} (Ai,...,A„). (2.5) 

If we solve the last equation with respect to the Lagrange multipliers we can express them in terms of the measured 
meanvalues 

K = K{G 1 ,...,G n ). (2.6) 



Whe n we substitute the Lagrange multipliers (2.6) into the expression for the generalized canonical density operator 
( |2.3|) we obtain the explicit expression for the reconstructed (estimated) density operator. 

As a typical example of the application of the MaxEnt principle we can consider a measurement of a single-mode 
electromagnetic field, modeled as a harmonic oscillator. Imagine, that as a result of the measurement we know the 
mean photon number n in the given field mode. Certainly, there are (infinitely) many quantum states of a single-mode 
electromagnetic field (e.g., a Fock state, a coherent state, a squeezed state, etc.) with the given mean photon number. 
So the question is: Which is the best (most reliable) estimation of the measured state given the mean photon number 
is known? The mean photon number is in some sense the least available information about the measured state. 
Consequently, the state is least determined. On the other, pure states are completely determined which is reflected 
by the fact that they have a zero von Neumann entropy. Therefore we expect that the most reliable reconstruction 
in the given case is a statistical mixture which is determined just by a single parameter - the mean photon number. 
It is well known that a statistical mixture which is parametrized just by a single parameter is a thermal state. This 
statistical mixture of Fock states is characterized by a temperature, or the corresponding mean photon number. In 
addition, for a given temperature (mean photon number) the thermal state exhibits the largest von Neumann entropy. 
Consequently, from the Jaynes principle of the maximum entropy it follows that if from a measurement only a mean 
photon number is known that the most reliable estimation of the measured state is the thermal state. 

The MaxEnt principle is not the only criterion how to choose an appropriate density operator among those 
which fulfill the constraints (p.l|). Based on an intuition or some addition a priori knowledge one can apply other 
criteria. For example, the maximum likelihood principle has been adopted successfully for estimation of quantum 
states p"2fl . Although this reconstruction scheme can result in nonphysical estimations (e.g., density operators which 
are not properly normalized, etc). In general, a consistent reconstruction scheme has to avoid nonphysical results 
(e.g. occurrence of negative probabiliti es). In the MaxEnt reconstruction a physical estimate is guaranteed by the 
canonical form of the density operator (2.3). The MaxEnt principle is the most conservative assignment in the sense 



that it does not permit one to draw any conclusions not warranted by the data. From this point of view the MaxEnt 
principle has a very close relation (or can be understood as the generalization) of the Laplace's principle of indifference 
which states that where nothing is known one should choose a constant valued function to reflect this ignorance. Then 
it is just a question how to quantify a degree of this ignorance. If we choose an entropy to quantify the ignorance, 
then the relation between the Laplace's indifference principle and the Jaynes principle of the Maximum Entropy is 
transparent, i.e. for a constant-valued probability distribution the entropy takes its maximum value. 

The MaxEnt reconstruction has been applied for various quantum systems, such as light field mode, spin systems 
jL]J. In what follows we adopt it for the reconstruction of vibrational states of neutral atoms. We assume the 
experimental setup realized by the group of C. Salomon at the ENS in Paris Qg. 
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III. RECONSTRUCTION OF MOTIONAL STATES OF NEUTRAL ATOMS 



Recently, experimental manipulations of motional quantum states of neutral atoms have been reported by the group 
of C. Salomon in Paris |||||. Cold Cs atoms can be cooled into specific quantum states of a far detuned ID optical 
lattice. The optical lattice is induced by the interference of two laser beams. Along the vertical z axis a periodic 
potential of "harmonic" micro-wells is produced with a period of 665 nm and with an amplitude of about 0.2 M. 
The vertical oscillation frequency in a micro-well at the center of the trap is w z /27r = 85 kHz. The corresponding 
ground state has the rms size Azo = y / h/2muj z « 21 nm and Ap /m = yjfiuj z /2m « 11 mm/s is its rms velocity 
width. The trapped cloud of neutral Cs atoms has a nearly Gaussian shape with a vertical rms size A£ = 53 p,m. 
With the help of deterministic manipulations the neutral atoms can be prepared in non-classical ID motional states 
along the vertical axis such as squeezed states, number states, or specific superpositions of number states j^j. The 
measurement of the prepared quantum state p is performed as follows: The system is evolved within the harmonic 
potential during the time r. Then the lasers are turned off and the system undergoes the ballistic expansion (BE). 
After the time of flight T — 8.7 ms a 2D absorption image of the cloud is taken in 50 (is with a horizontal beam ||. 
Integration of 2D absorption images in the horizontal direction gives us the spatial distribution along the vertical z 
axis. Therefore we will consider only ID quantum-mechanical system along the vertical axis. 

To confirm that a desired quantum state has been obtained (engineered) one can compare the spatial distributions 
along the vertical axis with the predicted ones. The coincidence of these spatial distributions is a necessary but not 
the sufficient requirement. A complete verification of the fidelity of the preparation of desired quantum states requires 
a quantum state reconstruction procedure. In order to perform this task we adopt the MaxEnt principle |jTl| . To do 
so we utilize a close analogy between quantum homodyne tomography jlij ] and the BE absorption imaging for the 
case of the point-like cloud (with the rms size equal to zero). 



A. Quantum tomography via MaxEnt principle 

Quantum tomography is based on the inverse Radon transformation of the the measured probability density distri- 
butions Wp(xg) for rotated quadratures Xg = (ae~ ie + a^e* 6 *) |l3] , |l4j| . These distributions can be represented as a 
result of the measurement of the continuous set of projectors |xe)(a;e|. Based on the measurement of the distributions 
Wp(xe) for all values of 9 £ [0, it] we can formally reconstruct the density operator according to the formula Jl5|] 
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d6 I dx e \xg)(xg\X(xg) 

J — oo 

where the Lagrange multipliers X(xg) are given by an infinite set of equations 



(3.1) 



Wp(xg) = \f2~n{xg\p r \xg)\ Vxe e (-oo, oo). (3.2) 

If the distributions Wp(xg) are measured for all values of xg and all angles 9 then the density operator p r is reconstructed 
precisely and is equal to density operator obtained with the help of the inverse Radon transformation or with the help 
of the pattern functions (for more details see |i~5|). 

In practical experimental situation (e.g., see the experiments by Raymer et al. [|l6| and by Mlynek et al. [pT[) 
it is impossible to measure the distributions Wp(xg) for all values of xg and all angles 9. What is measured are 
distributions (histograms) for finite number Ng quadrature angles 9 and the finite number N x of "bins" for quadrature 
operators. This means that practical experiments are associated with an observation level specified by a finite number 
of observables 

Fjk = \xS)(x^\ (3.3) 

with the number of quadrature angles equal to Ng and the number of bins for each quadrature equal to N x . Wc 
can therefore assume that from the measurement of the observables Fjk the meanvalues Fjk are determined (these 
meanvalues correspond to "discretized" quadrature distributions). In addition it is usually the case that the mean 
excitation number of the state is known (measured) as well. 

The operators Fjk together with h form a specific observation level corresponding to the incomplete tomographic 
measurement. In this case we can express the generalized canonical density operator in the form 



|exp 



N m Ng 

■^n-EEW I ;:U) 



j=l k=l 
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The knowledge of the mean photon number is essential for the MaxEnt reconstruction because it formally regularizes 
the MaxEnt reconstruction scheme (the generalized partition function is finite in this case). 



B. Motional states of atoms via MaxEnt principle: Formalism 



In the quantum homodyne tomography the probability distributions are measured for the rotated quadrature 

are related to the position 



a and 



it 



operators x$. The annihilation and creation operators of motional quanta, 

and momentum operators, z and p, via expressions z — + and p — — respectively. The angle 

9 of the quadrature operator corresponds to uj z t and vertical "cuts" of the absorption images (taken after the BE) 
can be associated with quadrature probability distributions. However, for a real physical situation with a nonzero 
rms size of the cloud the vertical "cuts" of absorption images correspond to a coarse-grained quadrature probability 
distributions. In particular, the vertical cuts of measured absorption images (taken in 2D) give us (after integration 
along the horizontal direction) the spatial distribution along the vertical axis. The spatial distribution along the 
vertical z axis can be expressed as 



F T (z)=T~ 1 J Fo(^)Pr((z-^)/T)d^, (3.5) 

where -Fb(£o) is the initial spatial distribution of the cloud in the z-direction (i.e., a Gaussian distribution with the 
rms size A£o)- The function P T {v) denotes the velocity probability distribution of the measured quantum state which 
has been evolved for time r in the harmonic potential before the BE, i.e. 



P T (v) = \( V \^(T))f 



\i>(r)) = tT(T)|^(0)>. 



(3.6) 



Here U(r) = exp(— iiTr/Ti) represents the time-evolution operator for the harmonic oscillator with the Hamiltonian 
H = p 2 /2m + mui 2 z 2 /2. Now we can treat the measured "cuts" as the mean values of specific observables: F T (z) = 
Tr [pF T (z)\. In practice just a few discrete times tj (J = 1, . . . , N T ) are considered and the z coordinate is discretized 
into the bins Zk (k = —N Z) . . . , N z ) of a given resolution Az. The set of operators which enters the equation ( |3.4| ) for 
the MaxEnt reconstruction then takes the form 



F Tj (z k ) = T- 1 / Fo^uHtj) 



Zk — Co \ / Zk — £o 



T 



T 



U(tMo (j = l,...,N T ;k = -N z 



,N Z ). 



(3.7) 



We have already commented that the operator of mean phonon number n is added to the set of observables {F Tj (zu)}- 
Knowledge of the mean excitation number n is essential in the case of an incomplete set of observables . Knowledge 
of the mean excitation number leads to a natural "truncation" of the Hilbert space. The inclusion of the mean phonon 
number into the MaxEnt reconstruction scheme does not represent its limitation as the mean energy represents one 
of basic characteristics of any system which should be inferred from the measurement. 

The experimental "cuts" of the BE absorption images [F Tj (z)] can be taken at few selected times, for example 
u) z Tj — 0, 7r/4, 7r/2, 37r/4 (N t = 4). To perform the reconstruction we have to determine the Lagrange multipliers 
{A^fc} an d A„ associated with {F Tj (zk)} and ft, respectively, in the expression for the generalized canonical density 
operator ( ft.4| ). The Lagrange multipliers can be determined via the minimization of a deviation function AF with 
respect to the measured data, i.e. 



AF = Y^ w jlk [F T] {z k ) - Tr [p r F Tj (z k j) }' 



w n {n - Tr (p r n)} 



(3.8) 



Here {wj^} and w fi represent positive weight factors for particular observables. Without any prior knowledge about 
the state we can take for simplicity Wij = 1. The weight factor w ft associated with the mean phonon number can 
be chosen according to our preference either to fit better the "cuts" of the BE images or the mean phonon number. 
In the case of the perfect measurement and the complete reconstruction the result has to be independent of the 
choice of the weight factors (in this case we can take Wn — 1). The weight factors could be also associated with the 
prior information about the dispersion of the measured observables. In particular, the weight factors can be taken as 
w v ~ <j~ 2 to reflect the knowledge of variances a v for the measured observables G v . When the mean values of the 
observables for the MaxEnt estimate p r fit within desired interval G v ± a v then contributions of the observables to 
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the deviation function AF are of the same order (~ 1). However, in our case we do not assume the knowledge of 
variances for the measured discretized probability distributions (taking w„ = 1). 

Once the Lagrange multipliers are numerically fitted, the result of the reconstruction - the generalized canonical 
density operator p r - can be visualized, for example, via the corresponding Wigner function fl8|| which can be defined 
as a particular Fourier transform of the density operator p of a harmonic oscillator expressed in the basis of the 
eigenvectors \q) of the position operator q: 

/oo 
dC(q-Cm P \q + C/2)e^. (3.9) 
-oo 



C. Numerical simulation 

To test our reconstruction procedure let us consider the reconstruction of the Wigner function of the motional 
quantum state IV'(O)) = (|0) + |l))/\/2 of Cs atoms trapped in the optical lattice. This kind of states has been 
demonstrated in recent experiments ||. We assume the following setup parameters: uj z /2tt = 80 kHz, the rms size 
of the ground state Az = 22 nm, the rms velocity width Apo/m = 11 mm/s and the rms width of the cloud of the 
atoms about 60 pm. Before BE (with BE time T = 8.7 ms) the atoms evolve within the harmonic trapping potential 
for t = 0, 1.6, 3.2, 4.8 ps. As the input for the reconstruction via the MaxEnt principle four vertical "ideal" cuts of the 
BE absorption images are taken as shown in insets of Fig. 1 . In addition, for the phonon number operator n which is 
included in the set of measured observables (see discussion above) we assume the mean value h = 0.5. The result of 
the ideal reconstruction is shown in Fig. 1. The fidelity of the measured and the reconstructed quantum states is close 
to unity which means a perfect reconstruction with AF = 10~ 10 , entropy S = 10 -7 , Ap = 10~ 8 has been achieved. 
Here Ap = n \ (p — p r )mn\ 2 denotes a deviation of the original and reconstructed density operators. 

r 2 =0 ^/s r 2 =1.6 [is 




FIG. 1. Numerical simulation of the reconstruction of the Wigner function of the motional quantum state (|0) + |l))/v / 2 of 
Cs atoms trapped in the optical lattice (assuming lu z /2ti — 80 kHz, the rms size of the ground state Azo = 22 nm and the rms 
velocity width Apo/m — 11 mm/s). As shown in insets for the reconstruction via the MaxEnt principle four vertical cuts of 
the absorption images (with BE time T = 8.7 ms) have been taken. The histograms in the insets correspond to the measured 
data while the solid lines are obtained from the reconstructed Wigner function (i.e. they correspond to reconstructed marginal 
distrubutions) . Before BE the atoms evolve within the trapping potential for the times r = 0, 1.6, 3.2, 4.8 ps. In addition, the 
mean number of motional quanta n = 0.5 and the rms width of the cloud of the atoms about 60 pm have been assumed. 

Obviously, in a real measurement the measured values are always fluctuating around the exact ones due to an 
experimental noise. Therefore we simulate a non-ideal measurement introducing random fluctuations to the measured 
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values of observables. It means that instead of the ideal values F Tj (z k ) we use for the MaxEnt reconstruction procedure 
the fluctuating ("noisy") values 

T' T . (z k ) = F Tj (z k ) + ^ (F T] (z k )) 1/2 . (3.10) 

Here ?; is a relative-error parameter which characterizes the quality of the measurement and {£j,/c} represents a 
Gaussian noise for observables. The result of the reconstruction is shown in Fig. 2 for rj = 0.1. Noisy mean values of 
the observables are shown in the insets. Despite of a significant relative error the reconstruction is almost perfect with 
the fidelity of the measured and the reconstructed states still close to one (AF — 0.16, entropy S — 0.01, Ap = 0.05). 
The minimum value of the deviation function AF = 0.16 can serve also as a measure of the imperfection of the given 
measurement (due to a technical noise) |2(J . 

r 2 =0 r 2 =1.6 




FIG. 2. Numerical simulation of the reconstruction of the Wigner function of the atomic motional quantum state 
(|0> + l})/\/2 for the same settings as in Fig. 1. Four vertical cuts of the absorption images taken for reconstruction (see 
the insets) are fluctuating randomly around their ideal values shown in Fig. 1 with the relative error r\ — 0.1. In addition, the 
mean phonon number fi — 0.6 has been considered. 

A typical non-classical state which we can utilize for a further test is the even coherent state M e ( | a) + 1 — a) ) which is 
a superposition of two coherent states with opposite phases For the amplitude a = \/2 we obtained AF = 10~ 8 , 
the entropy S — 0.026 and Ap = 10~ 4 (under assumption that the exact mean phonon number n — 1.928 is known). 
In the case of the imperfect measurement with rj = 0.1 the reconstruction leads to AF = 0.14, entropy S — 0.13, and 
Ap = 0.06 for n = 2.09. The fidelity of the reconstructed and the measured states is in this case also close to one. 

In order to model a technical noise in the measurement we have been considered Gaussian fluctuations proportional 
to the square root of the mean values. It means that tails of the "cuts" of BE images do not introduce a significant 
error (compare insets in Fig. 1 and Fig. 2). However, in the current measurements the situation seems to be different 
and the fluctuations do not decrease with the amplitude of the expected values. 

The fundamental question in the context of the MaxEnt reconstruction of states from incomplete tomographic data 
is whether the quality of the reconstruction can be improved using additional data from subsequent time moments r 
and how many such time moments r are required for the complete reconstruction of the unknown state p. We have 
shown recently [fl5| that for the quantum tomography just three quadrature distribution are sufficient for a complete 
reconstruction using the MaxEnt principle (in the case of the perfect measurement). This corresponds to the ideal 
case without the spatial dispersion of the cloud of atoms, i.e. the choice with u> z Tj = 0, tt/4, 7r/2 (N t = 3) is sufficient 
for A£o — ► 0. Obviously, in experiments with neutral atoms the spatial size of the atomic cloud is nonzero. However, 
in the case of the ideal measurement three BE absorption images associated with three "rotations" ui z Tj are still 
sufficient for a complete reconstruction of tested examples of quantum states. On the other hand, it seems that for 
higher mean phonon numbers the spatial distributions along the vertical axis which are directly determined from 
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absorption images should be known with improving precision (and on a wider interval of values as well) . In the above 
examples we have considered for convenience BE images for four "rotations" (N T = 4) which results in a very good 
reconstruction. 
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FIG. 3. Numerical simulation of the reconstruction of the Wigner function of the motional quantum state M e {\a) + \ — a)) 
with a — \[2 for the case of the ideal measurement. The mean number of motional quanta n = 1.928. Other settings are same 
as in Fig. 1. 



r 2 =0 /is 



r 2 =1.6 /is 



0.4 



0.2 



a 1 




0.4 

Jo.3 

•S 0.2 
I 

1 0.1 

.s 







-4 -2 2 4 
r 2 =3.2 /is 











0.4 




0.3 




0.2 




0.1 




0.0 





-400 -200 200 400 
z [/iin] 



-400 -200 200 400 
r 2 =4.8 fiS 




-400 -200 200 400 
z [/i.m] 





FIG. 4. Numerical simulation of the reconstruction of the Wigner function of the motional quantum state J\f e (\a) + | — a)) 
with a = \/2 for the case of noisy measurement with 77 = 0.1. The "measured" mean number of the motional quanta ft! = 2.09. 
Other settings are the same as in Fig. 1. 
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IV. RECONSTRUCTION FROM EXPERIMENTAL DATA 



In what follows we will apply the MaxEnt reconstruction scheme to the data obtained at the ENS in Paris 
Firstly we note that the unknown quantum state should belong to a Hilbert subspace which can be determined easily. 
Thus we can limit ourselves to the subspace spanned by Fock (number) states |0), . . . \N — 1). The upper bound 
on the accessible phonon number N is given by experimental limitations such as, for example, a feasible depth of 
micro- wells of the optical lattice and the validity of the harmonic potential approximation. For recent experiments N 
has been typically of the order of 10. This value is large enough to demonstrate the preparation of many non-classical 
states but on the other hand excludes highly squeezed states from a coherent processing. 

Let us consider the experimental arrangement used in Paris pl| with the parameters: uj z /2ir = 80 kHz, the rms 
size of the ground state Azq = 22 nm, the rms velocity width Ap$/m = 11 mm/s, the rms width of the cloud of the 
atoms about 60 /zm and BE time T — 8.7 ms. Initially the atoms are prepared in a well defined motional state \ipo) 
(e.g. in the vacuum state |0)). Then the optical lattice is switched off for the time period t\ during which the atoms 
evolve freely towards the state = exp(— itip 2 /2m)\ijjo)- Next, the optical lattice is again switched on for the time 
r during which the atoms evolve within the harmonic trapping potential. The measurement is performed after the 
BE time. The first two stages can be considered as the preparation of the state l^i)- After its "rotation" by lo z t 
(within the phase space of the harmonic oscillator) and the subsequent BE the absorption images are taken. 
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FIG. 5. The Wigner function reconstructed from the experimental data obtained at the ENS, Paris. The experimental 
setting is the same as for Fig. 1. From the experimental data we have inferred the mean number of motional quanta n « 1.0, 
while the reconstructed value is n w 1.1. Deviation of the measured and predicted values of observables is AF = 0.09 and 
entropy of the reconstructed mixture state is 5 = 1.0. 



The considered data are for the initial vacuum state which means that under ideal conditions a squeezed state 
= exp(— itip 2 /2m)\0) should be prepared. The vertical spatial distributions obtained from the measured 2D- 
absorption images are discretized into pixels (bins) with the pixel width 5.45 /im. The optical density of each pixel 
is averaged in the horizontal direction in which the absorption images are divided into 50 rows, each 3.9 p,m wide 
(these rows cover the size of the cloud in the horizontal direction). For the reconstruction via the MaxEnt principle 
four vertical spatial distributions for "rotation" times r = 0, 1.6, 3.2 and 4.8 /is are taken. The selected times roughly 
correspond to rotations within the phase space by lu z t — 0,7r/4, 7r/2 and 37r/4, respectively. Unfortunately, the 
mean excitation number n for measured state \tpi) was not measured explicitly in the experiment, therefore we have 
to estimate it as follows: During the free expansion period the rms size of the cloud increases by Aa; = poTi/m. 
The corresponding increase of the potential energy ^muil(Ax) 2 in units Tluj z gives us the increase of the number 
of excitation quanta with respect to the initial state IV'o)- F° r ti = A fj, and the initial vacuum it means n w 1. 
Experiments can be realized also for higher t\. For example, n = 8 fi leads to n ~ 4. However, as mentioned above, 
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such "squeezed" states with a significant contribution of higher phonon number states violate the underlying harmonic 
approximation for the potential. To keep a coherent control an an-harmonic part of the potential has to be taken into 
account. 

The result of the reconstruction via the MaxEnt principle is shown in Fig. 5. The deviation of the fitted and 
measured values is AF = 0.09 and the entropy of the reconstructed state S = 1.0. It means that the reconstructed 
state is a statistical mixture. We see a two peak structure which suggests that there is a mixture of two squeezed 
states coherently displaced from each other. It is caused by the fact that the vertical center of the cloud was not 
fixed in the experiment and it has to be determined by our fit for each measured BE absorption image separately. 
Assuming a priori knowledge that the Wigner function has a symmetric shape with respect to the origin of the phase 
space (i.e. there is no coherent amplitude) a Gaussian fit can be used to determine the center of the cloud for each 
vertical distribution. For states with a non-zero coherent amplitude the center of the cloud should be fixed already 
in the experiment. 

It turns out that the reconstruction results do not describe the squeezed vacuum state as was originally expected 
pl| . The main reason is that the mean phonon number was not measured directly in the experiment. It can be 
inferred only indirectly from the ideal case without any incoherence during preparation or measurement. As we 
discussed above, it is essential to include the information about the mean number of vibrational quanta into the 
MaxEnt reconstruction scheme. In optical tomography the analogous information about mean photon number can be 
obtained from distributions of two "orthogonal" quadratures. In our case it could correspond to two absorption images 
such that Lo z {rj — Tf.) = 7r/2. However, it would require a precise timing of the evolution within the harmonic trapping 
potential. Therefore the mean number of vibrational quanta should be determined in an independent measurement. 

Another problem arises from a slow convergence of anti-squeezed spatial distributions which are derived directly 
from the measured absorption images. In particular, the convergence of tails is too slow for those "rotations" which 
corresponds to anti-squeezed phases, i.e. r = 1.6, 4.8 /is (see insets of Fig. 5). The slow convergence is reflected by the 
presence of non-negligible backgrounds for Gaussian fits to these spatial distributions. If we eliminate (subtract) these 
backgrounds from the measured distributions the MaxEnt reconstruction gives almost the same Wigner function as 
in Fig. 5 but with a highly reduced deviation function AF = 0.02 (comparing to AF = 0.09 in Fig. 5). Such 
background in these absorption images can be caused by an incoherence associated with a violation of the the 
harmonic approximation. In fact, in our analysis we have neglected the change of the oscillation frequency along the 
z-axis. In recent experiments, the oscillation frequency decreases 10% from w 2 for micro-wells at the edge of the initial 
cloud. 



V. CONCLUSIONS 

Wc have presented a very efficient reconstruction scheme for the reconstruction of motional states of atoms based 
on the MaxEnt principle. The main advantage of the scheme is that it always results in reconstructions which are 
physical states (unlike in the case of the maximum likelihood estimation which can results in nonphysical estimations). 
Moreover, the scheme is very efficient in a sense that it requires just a small number of tomographic phase-space cuts. 

We have applied the scheme for a reconstruction of motional quantum states of neutral atoms. As an example have 
we analyzed the experimental data obtained by by the group of C. Salomon at the ENS in Paris and we reconstruct 
the Wigner function of motional quantum states of Cs atoms trapped in the optical lattice. In our analysis we have 
neglected the change of the oscillation frequency along z axis in recent experiments. The dispersion of the oscillation 
frequency is of the order of a few percent. This source of errors can significantly affect the quality of a quantum state 
preparation and its reconstruction. In addition, only up to first 10 bound states of micro-wells of the optical lattice 
can be approximated by a harmonic potential. It implies limits on coherent manipulations of quantum states. It 
means that states with a significant contribution of higher number (Fock) states cannot be prepared and manipulated 
in a controlled way. 
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